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We study the two-dimensional localization problem for (i) a classical diffusing particle advected by a 
quenched random mean-zero vorticity field, and (ii) a quantum particle in a quenched random mean- 
zero magnetic field. Through a combination of numerical and analytic techniques we argue that 
both systems have extended eigenstates at a special point in the spectrum, E c , where a sublattice 
decomposition obtains. In a neighborhood of this point, the Lyapunov exponents of the transfer- 
matrices acquire ratios characteristic of conformal invariance allowing an indirect determination of 
1/r for the typical spatial decay of eigenstates. 

PACS numbers: 46.10.+Z, 05.40. +j 
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In this paper we study two simple models for passive 
advection of a diffusing field: (I) the diffusion of a scalar 
density, n(x), advected by a quenched random velocity 
field, A(x), described by the Fokker-Planck equation pi: 



dm = C 



fpi 



DV 2 n - V • (An) 



(I) 



where D is the diffusivity; and (II) the random-flux 
model |^] for a non-interacting quantum particle prop- 
agating in a spatially-random, zero-mean magnetic field 
B = V x A = d x A y — d y A x , where A now denotes the 
vector potential, described by the Schrodinger equation 



-id t ip = Crrf = {(p - A) 2 + V}ip 



(II) 



ip being the (complex) quantum wave function, p = — iV 
the momentum operator, and V(x) the (scalar) potential 
disorder ||. 

Model (II) has received much attention recently in 
the context of the quantum Hall effect at filling factor 
v = I @' ^ n unres °l vec I question is whether the sys- 
tem has properties of a Fermi liquid, and in particu- 
lar, extended states. Previous work has addressed the 
energy-dependence of the localization length, £(E), mov- 
ing inward from the band edge, with authors arriving at 
opposite conclusions. The most careful numerical study 
of model (II) to date concludes that all states are local- 
ized [^J (b) , whereas others find a central band of extended 
states |£|(a). An analytic calculation using a replicated 
nonlinear sigma model with topological term [Q also ob- 
tains a band of extended states. 

Our purpose is two-fold. First, we explore the conse- 
quences of particle-hole symmetry at the band center of 
these models, E c , and describe numerical and analytical 
evidence for a divergent localization length at this point. 
Previous studies || |Q(6) did not allow for this symme- 
try at the band center |6[]. Second, we demonstrate that 
the properties of random flux that have drawn so much 
attention are exhibited by a much larger class of models, 
among them the passive scalar model (I). 



Magnetic field and vorticity are distinguished from po- 
tential fields by their transformation under time-reversal 
. Writing the velocity in model (I) as A = Vx + Vx(d, 
we observe that the curl-free part, Vx, like V, is even un- 
der time-reversal, whereas the divergence-free part V x <p 
(the source of vorticity to = VxA = — V 2 0), like B, 
is odd. A further physical similarity between the two 
models is that one expects transport to be dominated by 
the longest streamlines ; for the random-flux (passive 
scalar) model with vanishing mean magnetic field (vor- 
ticity), these rare streamlines run along the interfaces 
of opposing magnetic field (vorticity). The Laplacian 
enables fields to tunnel (diffuse) among distinct closed 
streamlines, and creates competition between advection 
by streamlines that can transport the field coherently 
over long distances, and diffusion that leads to destruc- 
tive interference |Io| . 

We study spatial decay of the eigenfunctions for lattice 
approximations to Cf p and C r f. Cf p is not self-adjoint, 
and its eigenvalues, z, occupy in general an area in the 
complex plane. For given z we use well-established nu- 
merical transfer-matrix methods and finite-size scaling 
to compute the localization length |l^] on a long 
strip. 

The real scalar field n is discretized on a square lattice 
of width L, length m. Various choices for boundary con- 
ditions in the transverse (L) direction will be discussed 
later, n and x are defined on nodes, A on links, and <f> on 
the nodes of the dual lattice. We define lattice difference 

i x and A 

where e 
tors. 



— V^x+e^ r * — • 

x, y, are orthogonal lattice basis vec- 



The velocity is then A = A x <fi + A + y, where 



A 



(-A 



y 



,A: 



represents the discrete curl. 
Defining the current J = DA + n — An, where [An]^ = 
ij4£[n x + n x+e J represents an average of the two ends 
of the link, the equation of motion is dtn + A~ • J = 0. 

Similarly, to discretize the random-flux Hamiltonian 
we define the lattice covariant derivative: 
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C rf = D+D- + V x 

The vector potential A has been defined on links in the 
same way as for the fluid; the scalar potential V* vanishes 
unless otherwise stated. 

Values of the field can be computed recursively using 
the 2L x 2L transfer-matrices, Wk, which yield values 
of ip or n in lattice column k + 1 and k given those 
in columns k and fc — 1. Lyapunov exponents are ex- 
tracted as logarithms of the 2L eigenvalues of the ma- 
trix (J^Mt^M)i/2m i n the limit m -> oo, where 
y\r( m ) = rj^^ Wk\ correlation lengths along the strip 
then correspond to their inverses. We define the scaled 
localization length by £,l{z) = 1/\l(z)L, where Xl(z) 
is the exponent smallest in magnitude. A critical or ex- 
tended phase occurs for those z where — > £00(2) 7^ 
for large L. 

We first describe results for model I for incompress- 
ible A = A~ x <fi, taking the </> x to be independent 
random variables distributed uniformly over an interval 
[— w, w]. When = the eigenvalues fill (uniformly, in 
d = 2) the real interval [-8D, 0]. For nonzero <p the den- 
sity of states broadens into a complex neighborhood of 
this interval. Fig. 1(a) displays using D = 1/4, 

w = 1, L = 32 and periodic boundary conditions. The 
peaks at z — 0, — 8D arise because the eigenfunctions, 
corresponding respectively to n.(x) = no and to uni- 
form antifcrromagnetic n(x) represent exact solutions for 
any width. The structure is symmetric about the line 
Re(z) = E c = — AD, a feature that originates in an ex- 
act particle-hole symmetry of Cf p - For L even, we divide 
the lattice into into its two equivalent antiferromagnetic 
sublattices, with n+ and n_ the restriction of n to the 
two sublattices, and A = ( ™ + ) ■ We can now express C f p 

in block form operating on A: 

where It, is the L x L unit matrix, Q = TQ T T~ 1 , T 
denotes transposition, and T time reversal. This trans- 
formation inverts the sign of the anti-symmetric parts 
of Cf p . The symmetry may be stated as follows: if 
A z is an eigenvector of Cf p with eigenvalue z, then 
A 2 = (T Z X* = (_"+» ) is an eigenvector with eigenvalue 

z = 2E C — z* . The same symmetry applies to the ran- 
dom flux operator, where E c — 4, z is real, and Q = QJ. 

At z = E c the two sublattices decouple, and the eigen- 
vectors of interest correspond to the zero eigenvalues of 
Q and Q. Furthermore, in the limit m — > oo, the Lya- 
punov exponents at E c must occur in degenerate pairs; 
this degeneracy is obtained numerically and disappears 
for any z / E c . The decoupling is also the source of the 
striking depression in £x, a t E c seen in Fig. 1 for both 
models. 



The depression and the degeneracy occur only for even 
L. For odd L, the boundary conditions mix the two sub- 
lattices. For periodic boundary conditions, £l(E) reaches 
(at E ~ —6.0 for the random flux model) a plateau of 
twice its degenerate value as one moves in from the band 
edges, and maintains that value at the band center. For 
odd L and free boundary conditions [as used in some nu- 
merical studies j^](a)], £l(E c ) diverges as m — > oo for any 
finite L. We obtain both degeneracy and divergence also 
for the "q" -models Q(a, c) where the fluxes are restricted 
to values 2im/q with q, n integers. These properties have 
not been identified before. 

Numerical evidence alone cannot distinguish an infi- 
nite localization length from a large but finite one. We 
now offer analytic arguments in support of a divergent 
correlation length at E c for the random-flux model with 
free boundary conditions. Let = ( n " fc ) represent the 
2L-component vector composed of the n x in columns k 
and k — 1. With an appropriate choice of gauge, we can 
write the transfer-matrix in the form Wk = ( ®* q L ) 

where Ofe is hermitian. Observe that WljWk = J where 
J = ( j l ) ; the set of matrices satisfying this iden- 
tity constitute a group. It follows that the eigenvalues of 
Wk occur in inverse conjugate pairs, /i, l//i*, as do the 
eigenvalues of any product of Wf-'s. 

The sublattice decomposition enables us to reor- 
ganize the components of fik in the form fik — 
col{n^, ^jt_i> n hi n k—i\- If i = 2Z + 1 is odd, the num- 
ber of components of will alternate with k between I 
and I + 1. The transfer-matrix now takes the block form 

Wk = (7 I) where w k = (*/ and w k = ( e j "'). 
This new form of the transfer-matrix only connects sites 
on the same sublattice and the 9k are no longer hermitian 
nor (for odd L) necessarily square. 

Because the ensemble of random matrices weights 8k 
and Oj. equally, we expect that the eigenvalues of the sub- 
matrix products ( w MW™))l/2m and ^(m)t !i (m)jl/2m 

are identical in the limit m — + oo, where w^ m > = IlfcLi w k 
and w^" 1 ) = IlfcLi^fc- It follows that the eigenvalues 
of the full transfer-matrix product occur in degenerate 
pairs; this degeneracy is observed numerically. Since the 
full transfer-matrix product is a group element, we de- 
duce that for odd L, there must be a pair of eigenvalues 
with modulus unity, one from each of the two submatrix 
products, yielding a divergent £l(E c ). 

For even L, an eigenvalue of modulus unity is not ex- 
pected for finite L, and we instead argue that a pair of 
eigenstates exists at E c in the thermodynamic limit. To 
make further progress we turn our attention from the 
transfer-matrix W to the operator C r f itself, and exploit 
its special form (4) at E c : the singular value decompo- 
sition of Q. (For convenience, we translate the of en- 
ergy to E c in this discussion). For a lattice with an odd 
number of sites, TV, Q is not square so that it has a non- 
trivial kernel and is a two-fold degenerate eigenvalue 
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of C r f. As a consequence of the singular value decom- 
position, adding a new lattice site (reversing the parity 
of N) can never increase the magnitude of the smallest 
non-zero eigenvalue. Because the randomness in Q can 
be expected to remove any accidental exact or near de- 
generacy, we anticipate rather that the magnitude of the 
smallest eigenvalue above (and its particle-hole conju- 
gate below 0) diminishes as N — ► oo, yielding a degener- 
ate pair of eigenvalues. 

If we accept that in the limit N — > oo through even val- 
ues there are indeed two independent eigenfunctions at 0, 
they must take the form (± v ), u,v the left, right eigen- 
vectors of Q so that Q^u = Qv = 0. Because < u\v >^ 
in general, we see that the arbitrary relative phase of u 
and v implies a continuous 17(1) symmetry at E c . Such a 
continuous symmetry is known to play an important role 
in some closely-related random-matrix models. Wegner 
first observed the significance of the sublattice decompo- 
sition in a class of random-matrix models for localization 
in d = 2 |0. It was later noticed that the sublattice 
decomposition allows a new continuous symmetry, which 
contains in the n = replica limit a factor of £7(1). For 
lattice models with spin, this continuous symmetry leads 
to a divergent DOS and localization length |L3| |. Our 
numerics indicate a finite DOS at E c , a result that in 
general is not inconsistent with a divergent correlation 
length. We are pursuing an analogous replicated field 
theory calculation for our model; results so far are con- 
sistent with the existence of this symmetry jl4| . 

We turn now to the neighborhood of the band center. 
Assuming finite-size scaling, we expect £l(z) to be inde- 
pendent of L (though not of z: see [la] ) as L — > oo in 
the regime of extended states. In Fig. 2, we show £l{z) 
for several values of z and L. Numerical values for £l(z) 
are more or less independent of L when z is on the lines 
Im(z) = or Re(z) = E c . 

For E c and its neighborhood we have examined the 
entire Lyapunov spectrum, and find that after appropri- 
ate scaling (see Fig. 3) the spectra for distinct L collapse 
onto a single curve. In addition, for small n/L, ratios of 
Lyapunov exponents take the form 

X n JXl = 2n - 1 (5) 

where \ 2 L denotes the j-th largest positive Lyapunov ex- 
ponent. This relation is remarkably universal: it con- 
tinues to obtain for both models (I) and (II) even when 
continued to complex A (this allows an approximate in- 
terpolation between the two models), and also when the 
B x (rather than the A£) in model II are taken to be 
independent random variables. As explained elsewhere 
|^5| , the form R n = (n + x)/x for the ratio of Lyapunov 
exponents suggests that the single-particle Green's func- 
tion is conformally invariant for a typical realization of 
the disorder, and typically decays as l/r v , where rj = 2x. 
Evidently, for our models r; ~ 1. 



So far we have discussed only divergence-free A. For 
model (I), the addition of a random \ leads to a real 
non-zero diagonal component of the discretized model at 
E c [as does scalar potential disorder V x for model (II) - 
note that x ma y always be removed by a gauge trans- 
formation]. The diagonal components of the disorder 
invalidate particle-hole symmetry and the sublattice de- 
composition. For small w, our numerics for these cases 
are nevertheless consistent with extended states, and the 
form x ~ 1/2 for the ratio still obtains. For a pure curl- 
free random velocity field (only \ nonzero), we find that 
all states are localized. 

We have done further numerical computations in order 
to isolate the property of the operators C f p and C r f re- 
sponsible for the apparent band of extended states @ . 
It is found that symmetric real random-matrix models 
with a sublattice decomposition display at E c degenerate 
extended states and x ~ 1/2, but away from E c all states 
are unambiguously localized, as is observed for the q = 2 
model [|)(a), and R n has no special form. The addition 
of a random, uncorrelated anti-symmetric matrix (real 
or complex) produces, in the neighborhood of E c , a band 
of states for which the Lyapunov ratios satisfy x a 1/2 
and £l has no detectable dependence on L. Recalling 
that the symplectic group constitutes the set of transfor- 
mations under which a bilinear anti-symmetric form is 
invariant, we conjecture that the usual assumption that 
the universality class of these models is described only 
by a unitary symmetry || M (b) may not be justified. 

We propose that chiral phase fluids, systems with large 
or divergent correlation lengths originating in a random 
anti-symmetric contribution to their dynamics, are a gen- 
eral phenomenon and share universal features. A full un- 
derstanding of these systems must answer the question, 
not so far resolved numerically, of whether the only ex- 
tended states are in fact at the band center, and if so, 
what sets the scale of the correlation length in the re- 
mainder of the band. Since our arguments are general 
in nature, we expect to see similar universal behavior in 
other systems, among them the kinematic dynamo pq| . 

ACKNOWLEDGMENTS 

We are grateful to E. Martinec, S. Abanov, L. 
Kadanoff, P. Wiegmann, D.P. Arovas, and S-C. Zhang 
for advice and discussions. We especially thank P. We- 
ichman for extensively revising our manuscript. This 
work was supported primarily by the MRSEC Program 
of the National Science Foundation under Award Num- 
ber DMR-9400379. Extensive use was made of the CRAY 
C98/8128 at SDSC. 



3 



[1] J-P. Bouchaud and A. Georges, Phys. Rep. 4, 127, 
(1990); J-P. Bouchaud, A. Comtet, A. Georges, and P. le 
Doussal, Europhys. Lett. 3, 653 (1987); M.B. Isichenko, 
Rev. Mod. Phys. 64, 961 (1992). 

[2] (a) D.Z. Liu, X.C. Xic, S. Das Sarma, and S-C. Zhang, 
Phys. Rev. B, in press (1995); V. Kalmeyer, D. Wei, D.P. 
Arovas, and S-C. Zhang, Phys. Rev. B 46, 3116 (1993) 
and references contained therein; (b) T. Sugiyama and 
N. Nagaosa, Phys. Rev. Lett. 70, 1980 (1993). 

[3] Remarks on some connection between these problems 
may be found in: B. L. Altshuler and L. B. Ioffe, Phys. 
Rev. Lett. 69, 2979 (1992); A. Comtet, A. Georges, and 
P. Le Doussal, Physics Letters B 208, 487 (1988). 

[4] S-.C. Zhang and D. P. Arovas, Phys. Rev. Lett. 72, 1886 
(1994). 

[5] A.G. Aronov, A.D. Mirlin, and P. Wolile, Phys. Rev. B 
49, 16609 (1994). 

[6] On the other hand, the 2-d random Dirac operator stud- 
ied by A.W.W. Ludwig, M.P.A. Fisher, R. Shankar, and 
G. Grinstein, Phys. Rev. B 50, 7526 (1994), exhibits a 
number of similarities to our models, including a particle- 
hole symmetric point, E c , with a pair of extended states. 

[7] For an explanation of time-reversal in these systems, see: 
R. Graham, Helv. Phys. Act. 59, 241 (1986). 

[8] D.S. Fisher, D. Friedan, Z. Qiu, S.J. Shenker, and S.H. 
Shenker, Phys. Rev. A 31, 3841 (1985); D.S. Fisher, 
Phys. Rev. A 30, 960 (1984). 

[9] S.A. Trugman, Phys. Rev. B 27, 7439 (1983). 
[10] J. Miller (manuscript in preparation) argues that the 
mapping between these two models, and the phase tran- 
sition itself, originate in a super symmetry of the models. 
[11] A. Crisanti, G. Paladin, and A. Vulpiani, Products of ran- 
dom matrices in statistical physics, (NY: Springer- Verlag, 
1993) and references contained therein. 
[12] D.J. Thouless, J. Phys. C 5, 77 (1972). 

R. Oppermann and F. Wegner, Z. Phys. B 34, 327 
(1979). 

[13] R. Gade, B398, 499 (1993); R. Gade and F. Wegner, 
Nucl. Phys. B360, 213 (1991); F.J. Wegner, Phys. Rev. 
B 19, 783 (1979). 

[14] J. Miller, unpublished. 

[15] see e.g. D.H. Lee and Ziqiang Wang, Transitions between 
Hall plateaus and the dimerization transition of a Hub- 
bard chain, preprint (1995). 

[16] Ya. B. Zel'dovich, S. A. Molchanov, A. A. Ruzmaikin 
and D. D. Sokolov, Sov. Phys. Usp. 30, 353 (1987) and 
references therein. 

FIG. 1. (a) Scaled localization length, £l(z), for model (I) 
and contour plot. Parameters are D = ±, w = 1, L = 32 
with periodic boundary conditions. The inset shows £l(z) for 
model (II) , with L = 32 and fluxes chosen independently and 
uniformly on [0, 2tv] (filled circles). Empty circles display the 
ratio Statistical error in the £l(z) is ~ 5%. 



FIG. 2. Scaled localization length, £,l(z)/S,16(z), plotted as 
functions of L for 

various z. +, X, square, * and circle represent respectively 
z = -1.0, -1.0 + 0.25i, -1.0 + 0.6i, -1.25, -1.25 + 0.25i. Error 
bars are displayed when they exceed the symbol width. 

FIG. 3. Comparison of the ratios, R n = A n /A„ , for vari- 
ous different L with equation (5) (dashed line). The largest 
correlation length corresponds to n = 0. A linear fit for small 
n/L yields x = 0.50 ± 0.02. To minimize the error in Rn we 
choose a small no for which the error is less than 1%. We 
collapse the ratios for different L by defining y = (n+l/2)/L 
and rescaling the ratio by (2no + l)/L, where no oc L. Param- 
eters are: z = -1.0 + 0.25i, D = \ and n = 8,4, 2 for L = 64 
(+), 32 (X) and 16 (*), respectively. For z = —1.0 a similar 
picture is obtained with x = 0.49 ± 0.02. The inset shows, 
for the same data set, the relative difference between the nu- 
merically obtained values and equation (5). The dashed line 
displays the statistical error. 
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Figure 1 




1 



Figure 1 inset 
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Figure 2 
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